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ABSTRACT 

We analyze the intrinsic polarization of two classical Be stars in the process of losing their circumstellar 
disks via a Be to normal B star transition originally reported by Wisniewski et al. During each of five 
polarimetric outbursts which interupt these disk-loss events, we find that the ratio of the polarization 
across the Balmer jump (BJ-I-/BJ-) versus the V-band polarization traces a distinct loop structure as a 
function of time. Since the polarization change across the Balmer jump is a tracer of the innermost disk 
density whereas the V-band polarization is a tracer of the total scattering mass of the disk, we suggest 
such correlated loop structures in Balmer jump-V band polarization diagrams (BJV diagragms) provide 
a unique diagnostic of the radial distribution of mass within Be disks. We use the 3-D Monte Carlo 
radiation transfer code HDUST to reproduce the observed clockwise loops simply by turning "on/off" 
the mass decretion from the disk. We speculate that counter-clockwise loop structures we observe in BJV 
diagrams might be caused by the mass decretion rate changing between subsequent "on/off" sequences. 
Applying this new diagnostic to a larger sample of Be disk systems will provide insight into the time- 
dependent nature of each system's stellar decretion rate. 

Subject headings: circumstellar matter — stars: individual (pi Aquarii) — stars: individual (60 Cygni) 



1. INTRODUCTION 

Classical Be stars are well known to be characterized by 
having gaseous circumstellar decretion disks which are fed 
from mass-loss from the ir rapidly rotating central stars 
(iPorter fc RiviniusI [2001 . As detailed in the review of 
Carciofi 1(120101) . the viscous decretion disk model (VDDM; 
Lee et all Il99ll ) is adept at explaining many of the ob- 



servational properties of Be star disks, including the ob- 
served Keplerian rotation of the disk, and is generally 
considered the most promising model to explain the Be 
phenomen on, although alternate model s also have been 
propo sed (jBiorkman &: Cassinellil Il993l : ICassinelli et al.l 
|2002|) . The mechanism(s) responsible for ejecting mate- 
rial from the stellar photosphere in to a disk is also un- 
kown, a lthough both observational (iRivinius et al.l 1998t 
Neiner e t al. 200^) and theoretical ([Ando 19861 ICranmeij 



20091) studies suggest that non-radial pulsations might act 



to feed at least some of these decretion disks. 

While some classical Be stars exhibi t observ a tional 
evidence of la rge-scale asymmetries (lOkazakil Il997t 
IVakih et all [1998.: Stefi et al.., 20091) interpreted as arising 
from one-armed density waves " (IOkazakii[l99l . the gen- 
eral gas disk density has typically been m odeled by a very 
simple axisy mmetric power law (see e.g. lBiorkman[ll997t 
lPorteil[l99l . Rec ent contemporaneous op tical-infrared 
(IR) spectroscop ic (IW isniewski et al.ll2007a^ ■ polarimet- 
ric (ICarciofi et a l. 2007 ), a.nd optical -IR interferometric 
studi es (iTvclIeret al.ii2006t iGies et al.|[2007t iTvcner et all 
[20081 iPott et al.l 120101) ■ which are each most sensitive to 



different physical regions of disk, demonstrate how the 
radial-dependence of the gas density in these disks can 
be observationally constrained. Such works offer promise 
for testing the appropriateness of the single power-law 
adopted by most models. 

Some classical Be stars c an exhibit stable decretion 
disks for decades (e.g C Tau; IStefl et al]|2009t ): however, 
they are also known to experience aperiodic "Be to nor- 
mal B to Be" transitions whereby they lose (and subse- 
quently reg enerate) all observational signatures of hav- 
ing a disk (lUnderhill fc DoazanI [T983 IClark et al.ll2003l: 
IVinicius et al.ll2006[ ). The frequency of these transitions 
is not well constrained by observations, although the dis- 
covery of 12 new transient Be stars in multi-epoch study 
of eight open clusters suggests th at these major events 
are not rare (jMcSwain et al.l [20091 ). These systems rep- 
resent ideal testbeds to diagnose the fundamental mech- 
anism which drives disk formation in Be stars precisely 
because they are known to be actively losing (or gaining) 
a disk. Studying such disks with techniques capable of 
diagnosing the radial-dependence of the gas density, espe- 
cially in the inner-most regions of the disk, could enable 
an enhanced understanding of how these disks form. 

In paper I of this series (jWisniewski et all |2010| ) , we 
analyzed ~15 years of spectropolarimetric observations 
of the classical Be stars 60 Cygni and tt Aquarii which 
covered one disk-loss episode in each star, and discussed 
the time-scale and overall evolution of these events. In 
this paper, we present a first look at the behavior of 



^ Department of Astronomy, University of Washington, Box 351580 Seattle, WA 98195, USA, jwisnie@u. washington.edu, zhd@uw.edu 
^ NSF Astronomy & Astrophysics Postdoctoral Fellow 

^ Ritter Observatory, Department of Physics & Astronomy, Mail Stop 113, University of Toledo, Toledo, OH 43606, 
karen.bjorkman@utoledo.edu, jon@physics.utoledo.edu 

* Instituto de Astronomia, Geofi'sica e Ciencias Atmosfericas, Universidade de Sao Paulo, Rua do Matao 1226, Cidade Universitaria, 05508-900, 
Sao Paulo, SP, BRAZIL, carciofi@usp.br, xhaubois@astro.iag.usp.br 

® Space Astronomy Lab, University of Wisconsin-Madison, 1150 University Avenue, Madison, WI 53706, meade@astro.wisc.edu 
^ Faculty of Engineering, Hokkai-Gakuen University, Toyohira-ku, Sapporo 062-8605, Japan, okazaki@elsa.hokkai-s-u.ac.jp 



2 



the intrinsic polarization during these events and detail a 
unique new diagnostic which effectively traces the gas den- 
sity in the inner-most region of the disk. We investigate 
these "polarization-loop" signatures with a well vetted 3-D 
Monte Carlo NLTE code and present representative model 
runs which reproduce the observations and support our in- 
terpretation of this phenomenon. Finally, we outline our 
plans to perform more detailed modeling of these signa- 
tures. 

2. OBSERVATIONS 

TT Aquarii was observed 127 times between 1989 August 
8 and 2004 October 10 and 60 Cygni was observed 35 times 
between 1992 August 3 and 2004 September 26 with the 
University of Wisconsin's (UW) HPOL spcctropolarimc- 
ter, mounted on UWs 0.9m Pine Bluff Observatory (PBO) 
telescope. The full details regarding the observation, re- 
duction, and calibration of these data, along with basic 
properties of the total (interstel lar -I- intrinsic ) polariza- 
tion, are presented in IWisnicwski et al.l (|2010[ ). hereafter 
referred to as paper I or WIS2010. In this paper we an- 
alyze the intrinsic polarization of these data, which were 
obtained by removing the interstellar polarization com- 
ponent describ ed in Table 5 of WIS2010 using a modifie d 
Serkowski law (jSerkowski et al.lll975l:IWilking et al.lll982D . 
We note that the residual instrumental systematic errors 
depends mildy on the date of the observations and ranges 
from 0.027-0.095% in the U-band, 0.005-0.020% in the V- 
band, and 0.007-0.022% in the I-band. 

3. RESULTS: A NEW DISK DENSITY DIAGNOSTIC 

Spectropolarimetry can provide insight into t he density 
struc ture of circumstellar envelopes (see e.g. iBiorkmag 
l2000t ). Pre- or post-scattering absorption of these 
photons by hydrogen in the disk can imprint the 
wavelength-dependent opacity signature of hydrogen in 
the p olarization if the density of absorbers is signifi- 
cant ([Wood et al.l Il996"bl) . producing a "sa w-tooth" like 
wave l ength dependent polariz ation signature (jWood et al.l 
ll997HWisniewski et al.ll2007bD . Thus while the overall po- 
larization is a tracer of the electron number density, or ef- 
fective mass of a Be disk, the polarization change across 
the Balmer jump is a tracer of the largest densities of the 
disk, which are typically found in their innermost regions. 
With these basic principles in mind, we begin our analysis 
of the intrinsic spectropolarimetric dataset of WIS2010. 

WIS2010 detailed how the gradual disk-loss episodes of 
60 Cyg and tt Aqr proceeded over a timescalc of '^1000 and 
~2400 days respectively, and how these events were tem- 
porarily stalled by several polarimetric outbursts (left pan- 
els Figure [TJ sec also WIS2010). We examined the behav- 
ior of each of these events in standard Johnson broad-band 
filters, created by binning our spectropolarimetric data to 
reproduce the coverage of each filter, and in custom fil- 
ters such as one which probed the polarization across the 
Balmer jump. This Balmer jump (BJ-I-/BJ-) filter was 
created simply by ratioing data binned between the wave- 
length range 3650 - 4100 A (BJ-t-) and 3200 - 3650 A (BJ-). 

During polarimetric outbursts, we detect clear evidence 
that the polarization across the Balmer jump (BJ+/BJ-) 
exhibits a distinctive loop-like evolution when compared 
to the V-band polarization. We show examples of such 



loop-like structures in Balmer jump-V band polarization 
diagrams (hereafter BJV diagram) for both 60 Cyg (Fig- 
ure [TJ top middle panel) and tt Aqr (Figure [TJ bottom 
middle panel) . The outburst responsible for each of these 
loops are depicted in the top and bottom left panels of Fig- 
ure [TJ The duration of the depicted outburst in tt Aqr is 
^180 days and the evolution of the polarization across the 
Balmer jump traces a clockwise loop pattern as a function 
of V-band polarization. The duration of the depicted out- 
burst in 60 Cyg is ^770 days and while the polarization 
traces a clockwise loop pattern during the initial stages 
of the outburst, it traces a partial counter-clockwise loop 
pattern during the latter stages of the outburst. We find 
similar clockwise loop morphologies for one additional po- 
larimetric outburst in tt Aqr which lasted for ~160 days, 
and a counter-clockwise loop morphology for another tt 
Aqr outburst which lasted for ^^730 days. 

While these loops are very prominent in the BJV dia- 
gram (middle panels; Figure [TJ), we also explored whether 
these morphological features would be detectable if one 
simply had broad-band filters at their disposal. We there- 
fore binned our spectropolarimetric data to reproduce the 
coverage of the Johnson U- and B-filters, and plotted 
our data on a (B-filter/U-filtcr) versus V- filter (hereafter 
BUV) diagram. As seen in the right panels of Figurc[TJ the 
same looplikc structures are visible in BUV diagrams (as 
BJV diagrams), albeit at lower amplitude. We therefore 
suggest that moderately high precision filter polarimetry 
could also provide the same effective diagnostic as afforded 
by our spectropolarimetric data. 

4. DISCUSSION 
4.1. Modeling the Observed Phenomenon 

As mentioned in the introduction, the viscous decre- 
tion disk model (VDDM) is currently the most promising 
candidate to explain the structure, formation and evo- 
lution of Be disks. While most theoretical approaches 
so far considered the case of a constan t mass decretion 
rate in the quasi-steady sta te limit (e.g.. lBiorkman|[T997l : 
iBiorkman fc Carciofil l2005l ). more recent studies investi- 
gated the temporal evolution of the disk surf ace den- 
sity fed by an arbitrar y mass decretion history (jOkazakil 
[20071 1 Jones erall[2008h . Here we qualitatively investigate 
whether the VDDM can account for the trends we see in 
BJV diagrams. 

In order to theoretically reproduce the observed trends 
shown in Figure 1 wit h the VDDM model we used a 1-D 
hydrodynamical code ()Okazakil2007l ) to compute the time- 
dependent surface density of the disk. This code solves 
the viscous diffusion problem given a prescription for the 
stellar mass decretion rate, and a val ue for the disk kine- 
matic viscosity (the a parameter of IShakura fc SunvaevI 
|1973() . The surface density for chosen epochs of the disk 
evolution is then fed to our radiative transfer code hdust 
(jCarciofi &: Biorkmanll20d6t ) that is capable of turning the 
structural information thus provided into astrophysical ob- 
servables, such as emergent spectrum or intensity maps on 
the sky. 

The correlation loops observed in the BJV diagram can 
be qualitatively reproduced assuming a prescription for 
the mass decretion rate that involves alternating cycles of 
activity (mass decretion on) and quiescence (mass decre- 
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tion off). One such an example is shown Fig. [21 Start- 
ing from no preexisting disk, this sample model assumes 
a 3-year long period of activity followed by a 3-year long 
quiescence. This 6-year cycle were repeated many times. 
In Fig. [5] we show results covering the period between 32 
and 39 years after the beginning of disk formation. The 
polarization forms a clockwise loop in the BJV diagram 
that can be described as follows. At the end of the active 
phase, the star had built a large and dense circumstellar 
disk (phase 1). When the mass decretion stops, the inner 
disk quickly reaccretes back onto the star; this causes a 
fast drop of B J size and the curve follows a track towards 
the bottom-left of the BJV diagram (phase 2). What fol- 
lows is a slow secular dissipation of the entire disk along 
which the BJ size changes little (the inner disk having 
already reached very low densities) but the V-band polar- 
ization diminishes as the disk mass decreases (phase 2 to 
3). When the next cycle of activity begins (phase 3) the 
inner disk quickly fills up again and the curve eventually 
reaches back the top of the BJV diagram. 

The detailed shape of the loop depends on several fac- 
tors: the viewing angle (Fig. the value of a, and the 
mass decretion history assumed. In the simple model 
shown here the loops nearly close, because the mass de- 
cretion rate assumed for each cycle is the same. The loop 
would not close if the mass decretion rate of subsequent 
cycles were to be different and/or if the length of the ac- 
tive/quiescent phases were irregular. We plan to systemat- 
ically explore a large range of mass decretion rate scenarios 
in a future publication (Haubois et al 2011, in prep). 



4.2. Comparison with CMD loops 

Ide Wit et al.l ()2006l ) studied the photometric variabil- 
ity of several hundred Be stars in the Magellanic Clouds 
and found ~100 stars whose photometric variations traced 
loop-like patterns in optical color magnitude diagrams 
(CMDs). Most (^^90%) objects traced clockwise loops in 
these CMDs while 10% traced counter-clockwise loops. 
Ide Wit et al.l (|2G06i ) suggested that clockwise loops were 



indicative of systems actively decreting material from their 
stellar surfaces whereas counter-clockwise loops were in- 
dicative of systems in which material was being re-accreted 
onto the central star. 

We note that our models of the loop-like structures we 
observe in BJV diagrams include the effects of both decre- 
tion of material from the central star and the re-accretion 
of material onto the central star. Hence the large-scale 
morphological differences we observe, i.e. clockwise ver- 
sus counter-clockwise loops in BJV diagrams, can not be 
simply attributed to differences in the radi al direction of 
motio n of material in the disk as invoked by Ide Wit et al.l 
(|2006t ). Rather, we speculate that the counter-clockwise 
loops we sometimes observe in BJV diagrams might be 
caused by a non-constant mass decretion rate or two 
nearby cycles which have significantly different mass de- 
cretion rates. Our future systematic exploration of the 
parameter space of our models (Haubois et al 2011, in 
prep) will enable us to test this speculative hypothesis, 
as well as other mass ejection scenarios, to explain these 
counter-clockwise loops. 



4.3. Implications and Future Applications of the 
Technique 

The stea dy-state (iBjorkman fc Carciofill2005l) and time 
dependent (jOkazakil 2007f ) surface densitv of Be disks is 
driven in large part by the ratio of the stellar mass- loss rate 
and the a parameter (both quantities setting the disk de- 
cretion rate), although recent modeling work has demon- 
strated that the disk temperature can also influence the 
surface density, especially in the inner disk re gions, when 
the stellar mass-loss rate is large (jCarciofi et al. 2008) . As 
ICarciofi et al.l (|2009[ ) noted, observationally constraining 
the stellar mass-loss rate is very challenging for systems 
whose disks are truncated by binary companions; more- 
over, constraining the detailed time depen dent decretion 
rate of non-steady state disk systems (e.g. iRivinius et al.l 
Il998f) is also challenging. The BJV (and broad-band BUV) 
diagnostic we have introduced in this paper offers one clear 
way to better diagnose the detailed time-dependent decre- 
tion rate of Be disk systems. Application of this tech- 
nique, when the requisite low-resolution blue optical spec- 
tropolarimetry (BJV diagram) or U- and B-band filter po- 
larimetry is available, to a larger sample of Be systems 
actively gaining/losing their disks would help elucidate 
the mechanism(s) responsible for triggering disk forma- 
tion in Be stars. We encourage the community to obtain 
this type of well time sampled polarimetry for the mid- 
2011 periastron passage of 6 Sco, as it would be provide 
a powerful diagnostic of the type of decretion outbursts 
which have been reported in prev ious periastron passages 
(jMiroshnichenko et al.ll200lll2003[) . Moreover, inclusion of 
simple V-band photometry in such analys is could help ob - 
servationally determine the a parameter (jCarciofi Il2010[ ). 

5. SUMMARY 

We have analyzed the intrinsic polarization of 60 Cyg 
and TT Aqr as they were in the process of losing their 
circumstellar disks via a Be to normal B star transition. 
During each of five polarimetric outbursts which interupt 
these disk-loss events, we find that the ratio of the po- 
larization across the Balmer jump (BJ+/BJ-) versus the 
V-band polarization traced distinct loop structures as a 
function of time, and suggest that this diagnostic provides 
unique insight into the radial distribution of mass within 
Be disks. We observed two clockwise and one counter- 
clockwise loop structures in Balmer jump-V band polar- 
ization diagrams (BJV diagragms) of n Aqr, while 60 Cyg 
exhibited a combined clockwise and counterclockwise loop 
(i.e. a "figure eight"). We use the 3-D Monte Carlo ra- 
diation transfer code HDUST to reproduce the observed 
clockwise loops simply by turning "on/off" the mass decre- 
tion from the disk, and speculate that counter-clockwise 
loop structures might be caused by the mass decretion rate 
changing between subsequent "on/off" sequences. Current 
and future explorat ion of the parameter space of our mod- 
els (|Haubois II2010I : Haubois et al 2011 in prep) will help 
efforts to identify the definitive origin of these counter- 
clowckwise loop structures. 
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Fig. 1. — The time evolution of the intrinsic polarization of 60 Cyg (top panels) and tt Aqr (bottom panels) is shown. 60 Cyg is characterized 
by a ~770 day polarimetric outburst denoted by large, open green circles (top left panel). A clock- wise loop pattern in 60 Cyg's Balmer 
jump-V band polarization diagram (BJV diagram; top middle panel) is observed during the early stages of the star's outbust, followed by 
a counter-clockwise loop during later stages of the outburst. This loop structure is discernible at lower amplitude when the polarization 
ratio of the Johnson B/U filters is considered as a function of the V-band polarization (top right panel). The time evolution of the intrinsic 
polarization of tt Aqr is characterized by several polarimetric outbursts, including a '^180 day event denoted by large, filled red circles (bottom 
left panel). A broad clock-wise loop pattern is seen in tt Aqr's BJV diagram (bottom middle panel); this loop structure is discernible at lower 
amplitude when the polarization ratio of the Johnson B/U filters is considered as a function of the V-band polarization (bottom right panel). 
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Fig. 2. — Upper panel: BJV diagram at two inclination angles (solid line, 70 degrees and dashed line, 90 degrees) for the mass decretion 
rate described in Sect. 4.1. The corresponding phase number is reported in both panels and observing epochs are indicated and counted in 
years. Lower panel: Temporal evolution of the V-band polarization and polarization across the Balmer jump for both inclination angles. The 
dotted line shows the mass decretion history arbitrarily scaled to the range of the graph. 



